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Wc report evidence of large, nonlocal correlations between two spatially separated normal metals 
in superconductor/normal- metal (SN) heterostructures, which manifest themselves a nonlocal volt- 
age generated in response to a driving current. Unlike prior experiments in SN heterostructures, 
the nonlocal correlations are mediated not by a superconductor, but by a proximity-coupled normal 
metal. The nonlocal correlations extend over relatively long length scales in comparison to the 
superconduncting case. At very low temperatures, we find a reduction in the nonlocal voltage for 
small applied currents that cannot be explained by the quasiclassical theory of superconductivity. 
We believe is a signature of new long-range quantum correlations in the system. 



PACS numbers: 74.45-fc, 03.67.Mn, 74.78.Na 

Electrons in spatially separated normal metals in mu- 
tual contact with a superconductor show correlations 
that are evidence of quantum entanglement due to their 
interaction with the Cooper pairs in the superconductor. 
[1-8] Experimentally, the correlations manifest them- 
selves as a nonlocal voltage that develops on one nor- 
mal metal in response to a current injected from another 
normal metal into the superconductor. In such a hybrid 
structure, three processes that contribute to the nonlocal 
voltage have been investigated. Two, crossed Andreev re- 
flection (CAR) and elastic co-tunneling (EC), arise from 
the interaction of quasiparticles in the spatially separated 
normal metals that arc mediated by the Cooper pairs 
in the superconductor, and decay spatially on the scale 
of the superconducting coherence length ^g, typically ~ 
100 nm for Al [4, 6]. The third process, charge imbalance 
[9, 10], is associated with the conversion of a quasiparticle 
current into a supercurrent in the superconductor, and 
decays over the charge imbalance length Ag. , which is of 
the order of a few microns [11-13]. None of these three 
nonlocal processes are expected in a normal metal with 
an induced superconducting proximity effect, where there 
is no superconducting order parameter. Consequently, it 
comes as a surprise that large nonlocal signals can in- 
deed be observed when a quasiparticle current is injected 
into a normal metal that is proximity coupled to two su- 
perconductors. The experimental manifestation of these 
nonlocal signals is very similar to that observed in NSN 
structures [11-18], although their origin appears to be 
quite different. 

The samples measured in this work had Au as the nor- 
mal metal and Al as the superconductor. The length L of 
the normal wire between the two superconducting films 
was designed such that the two superconductors were 
Josephson-coupled in the temperature range of interest. 
For diffusive conductors, this occurs at temperatures be- 
low Ec/kB, where = hD/LP' is the Thouless energy 
and D the electronic diffusion coefficient. The tempera- 
ture dependence of Josephson coupling is determined by 



the Thouless length Lt = \/hD/kBT, where T is the 
temperature. For this sample, which had D=110 cm^/s 
so that Ec = 11.7 AteV and Lt = 293 nm/\/T, a finite 
supercurrent was observed below ^0.7 K. Given the mul- 
tiprobe nature of the sample, we shall use the notation 
Rij,ki = dVki/dlij to denote the four-terminal differen- 
tial resistance, where the ac and dc current is applied 
between contacts i and j, and the resulting ac voltage 
measured between contacts k and I, with the contacts 
numbered as in Fig. la. 

Figure lb shows a measurement of the local differen- 
tial resistance i?i9,28 a-s a function of Idc at T—20 mK. 
The position of the peaks in this curve identify the crit- 
ical current Ic- The blue circles in Fig. Ic show the 
T dependence of I^, and the solid line is a fit to the 
functional form = BT^/"^ eK\){—A/ Lt) (where A and 
B are constants) expected for a SNS junction in the 
long junction limit [19] (A >> Ec, where A is the 
gap in the superconductor). As can be seen, the fit is 
quite good. Ic at base temperature is smaller than the 
value /cO = 10.82£'c/ei?7v predicted for a simple SNS 
junction[19] {Rn is the normal state resistance). For this 
sample, Ico&Rn / Ec ~ 0.56, with Rm = 4.56 Vl. However, 
in multiterminal structures Ic is suppressed [20] due to a 
modification of the induccid niinigap. If we assume that 
the measured /cO is related to an effective Thouless en- 
ergy E* by the same theoretical prediction for a simple 
SNS junction, we obtain E^ = 0.59 fieV. 

We now discuss the nonlocal mcasurcinients. Figure; 2a 
shows the nonlocal dV/dl at 20 mK for four different con- 
figurations, each with the current sourced through nor- 
mal electrode 3 and drained through a superconducting 
electrode 1. The overall shape of the resulting traces is 
remarkably similar to what is observed in previous mea- 
surements on NSN devices [16, 17]: at Idc = 0, dV/dl 
is finite and grows with current, resulting in a peak at a 
finite current of ~ 2.3 /xA, after which there is a sharp 
drop to negative values before it goes to zero at high bias. 
In the case of NSN structures, a similar peak followed by 
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FIG. 1: (a) False color scanning electron micrograph of the 
sample discussed in the text. Light areas are Au; the darker 
lines are Al. The numbers mark the contacts used for four- 
terminal differential resistance measurements. The size bar is 
500 nm. (b) Local differential resistance i?i9,28 as a function 
of the applied current Idc, showing the critical current of the 
SNS junction, (c) Blue circles show the measured critical cur- 
rent, determined by the position of the peak in the differential 
resistance in (b), as a function of temperature T. The solid 
line shows a fit to the expected temperature dependence for 
a long SNS junction of length L [19]. Open red circles are the 
finite values of Idc at which the nonlocal differential resistance 
J23i,49 (Fig. 2(a)) has its minimum. 



a dip was associated with Idc exceeding Ic of the super- 
conductor. Here, as discussed below, this feature results 
from part of Idc exceeding Ic of the SNS junction [16, 21]. 
dV/dl also decreases as the distance of the normal volt- 
age contact from the current injection node increases. 
While the overall shape of the traces is similar to that 
observed nonlocal NSN devices, there are some signifi- 
cant differences. In particular, there is a sharp dip in 
dV/dl near Idc = that is not present in nonlocal NSN 
measurements. 

In the nonlocal NSN experiments [16, 17], it was found 
that the zero bias differential resistance i?„;(0) and the 
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FIG. 2: (a) Nonlocal resistances as a function of dc current 
bias for 4 nonlocal configurations, -R3i,4g, -Rsi.sg, -Rai.eg and 
^31, 79- (b) Data of (a), with the curves for -R3i,49, J?3i,59, 
^31, 69 and -R3i,79 scaled so that their normalized peaks at 
±2.3 /iA match, (c) Nonlocal resistance at 2.3 ^A (black 
symbols) and at zero bias (red symbols), taken from (a), as a 
function of the distance between the respective voltage con- 
tact and the current injection point. The straight lines are 
linear fits to the data, (d) Nonlocal resistances -R3i,79, i?4i,79 
and -R5i,79, where the current injection terminal is changed, 
but the voltage contacts remain the same. 

peak at finite Idc decayed with the distance from the cur- 
rent injection electrode, but with different length scales: 
while the zero bias resistance decayed exponentially with 
^5 as expected from CAR/EC, the peak was associated 
with charge imbalance and found to decay linearly with 
Aq*. In the current experiments, ii„/(0) and the peak 
resistance also scale differently with distance. Figure 2b 
shows the curves of Fig. 2a scaled so that their peaks 
at finite bias match. With this scaling, the curves match 
over most of the range of current, except near zero bias. 
This shows clearly that i?„;(0) and the finite bias resis- 
tance scale differently with length, as was found for the 
NSN samples. In the NSN case, i?„/(0) decayed expo- 
nentially on a length scale of ^s: Here one might expect 
that -R„;(0) should decay exponentially with Lt- Figure 
2c shows that, contrary to expectation, both Rni{Q) and 
the peak resistance show an almost perfect linear depen- 
dence on the distance from the current injection point, 
although the slopes are different, reflecting the different 
scaling evident in Fig. 2b. 

What is the origin of this behavior? In NS structures 
with Josephson coupling between two superconducting 
electrodes, supercurrents and quasiparticle currents can 
coexist in a proximity-coupled normal metal over dis- 
tances much longer than ^5 [22-24]. As the supercon- 
ducting electrodes Si and S2 are Josephson-coupled at 
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low enough temperatures, they are at the same electro- 
chemical potential, which we take to be here for simplic- 
ity. A finite potential V applied to the current injection 
electrode A^i will drive a quasiparticle current Iqp into the 
proximity- coupled normal metal (Fig. 3a). Since Si and 
52 are both at zero potential, this quasiparticle current 
will split into two components: Iqpi will flow towards 5*1, 
and Iqp2 will flow towards ^2, the ratio Iqpi/Iqp2 being 
determined by the inverse ratio of the resistances of the 
normal sections between the current injection point and 
Si, 2- As ^2 is a voltage probe, no net current can flow 
into it. Hence Iqp2 must be balanced by a counterflow- 
ing supercurrent Is that flows from 5*2 to Si such that 
Is = —Iqp2 and Is + Iqpi = Iqp, in turn giving rise to a 
phase difference (j) between 5*1 and 5*2 . This situation will 
persist until Is exceeds Ic of the junction. Since Iqpi is 
always less than Iqp, the injected dc current Iqp at which 
this occurs is always greater than Ic- Evidence for this 
model can be seen by examining the value of current at 
which this occurs. The red open circles in Fig. Ic show 
the value of Idc as a function of T at which the minimum 
in the resistance at finite bias is observed in the nonlo- 
cal measurement (Fig. 2a). The ordinate axis has been 
scaled so that the data points lie on top of the measured 
Ic (blue circles) at high temperatures. At lower temper- 
atures, however, they show a weaker temperature depen- 
dence. This latter behavior has been observed previously 
by other groups [24] , and arises from the reduction of Ic 
due to the nonequilibrium quasiparticle distribution in- 
troduced by Iqp [22]. 

Iqp2 will result in a finite "nonlocal" voltage between 
the normal voltage contacts A^2,3 and the second super- 
conductor S2 that will be proportional the resistance of 
the normal wire between A^2,3 and ^2, plus the resis- 
tance of the interface NS2- Thus, one would expect to 
see the linear scaling shown in Fig. 2c. In addition, as 
the current injection point is moved closer to 6*2, one 
would expect to see an increase in the nonlocal voltage 
(and hence differential resistance) as the ratio Iqp2/Iqp in- 
creases. Figure 2d, which shows the nonlocal differential 
resistance for a fixed voltage lead but different current 
injection points, confirms this. 

Hence it appears that the current separation model 
describes well the nonlocal resistance that we observe. 
However, closer analysis of the current bias and tem- 
perature dependence of the nonlocal resistance reveals 
some significant discrepancies. Quantitative predictions 
for the differential resistance can be obtained by simulta- 
neously solving numerically the Usadel equations as well 
as the kinetic equations for the nonequilibrium distribu- 
tion functions [26] to obtain the nonlocal dVni/dlqp as a 
function of Iqp, and i?„;(0) as a function of temperature. 
The resulting curves are shown in Figs. 3b and 3c re- 
spectively. (Details of the calculations are given in the 
Supplementary Information.) There are many differences 
between the results of the simulations and our measure- 
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FIG. 3: (a) Schematic of the current separation model. The 
injected quasiparticle current Iqp splits into two currents, Iqpi 
and Iqp2 , which go towards the two superconducting contacts 
Si and S2 respectively. Iqp2 is compensated by a counterflow- 
ing supercurrent 7s that flows from 5*2 to Si. (b) Results of 
the numerical simulations based on the quasiclassical theory 
of superconductivity for the nonlocal resistances for the ge- 
ometry of (a) as a function of Iqp. The distance L' between 
the normal probe (on the normal metal) and Si is 0.4511/, 
0.578L, 0.696L and 0.8431/, where L is the length of normal 
metal between Si and S2. The position of A''i is fixed at 
0.196L from Si. These values are chosen to match the sam- 
ple geometry, (c) Calculated zero bias resistance as a function 
of temperature. Details of the simulation are described in the 
Supplementary Information. 



ments, but we shall focus here on the most significant 
discrepancy. Figure 3c shows that i?„;(0) is expected 
to increase as T decreases. This behavior is simply the 
well-known reentrance effect of a proximity-coupled nor- 
mal metal [27]: as the temperature is decreased below 
the transition temperature of the superconductor, the 
resistance of the normal metal first decreases, but then 
reaches a minimum at a certain temperature Tq, increas- 
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FIG. 4: (a) Temperature dependence of Rni{0) for all four 
nonlocal resistance configurations, (b) Nonlocal differential 
resistance for the closest configuration at a number of different 
temperatures, (c) Temperature dependence of — /?„;(0) for 
the configuration in (b) and temperature dependence of the 
local critical current Ic- (d) Data from (b) at 20 mK, 80 mK 
and 160 mK, plotted as a function of the nonlocal voltage V„t 
obtained by numerical integration. Dotted lines are guides to 
the eyes for the points where 20 mK and 80 mK data deviate 
significantly from 160 mK data. 

ing as the temperature is decreased further. For a single 
N wire connected to a S reservoir, Tq ^ ^Ec/ks [28]. For 
more complicated geometries, Tq may be modified, and 
a self-consistent calculation as we have done is required. 
These calculations show that experimentally we are in 
the regime T < Tq, with the resistance rising as the tem- 
perature is decreased, since Lt > L in the temperature 
regime of interest. The experimental data, however, show 
exactly the opposite trend. Figure 4a shows i?„;(0) for 
the 4 nonlocal configurations of Fig. 2a in the low tem- 
perature regime. The nonlocal resistance decreases with 
decreasing temperature, in direct contrast to the behav- 
ior expected from Fig. 3c. At higher temperatures, the 
nonlocal resistance drops sharply to zero, which is due to 
the vanishing of the supercurrent between the supercon- 
ducting electrodes, without which no nonlocal signal can 
be observed. 

A clue to the origin of the decrease in i?„;(0) with 
decreasing T can be seen in the temperature evolution 
of dV/dl vs. Idc, which is shown in Fig. 4b for one 
nonlocal configuration. Apart from the decrease in Ic 
with increasing T, the only major difference between the 
different temperature traces is the growth of the dip at 
zero bias. In order to study the temperature dependence 
of this feature, we plot —Rni{^) together with the mea- 
sured Ic of the SNS junction, with appropriate offset and 
scaling so that the data points coincide at higher tem- 
peratures. — i?„i(0) matches the exponential behavior of 



Ic at higher temperatures, but the two curves diverge at 
lower temperatures, with Ic saturating, but —Rni(0) still 
showing a strong exponential dependence. To investigate 
the energy scales involved, we integrate the nonlocal dif- 
ferential resistance to obtain dV/dl vs Vni- Figure 4c 
shows the 20, 80 and 160 mK from Fig. 4b plotted in this 
manner. The voltage at which the 20 mK curve deviates 
significantly from the 80 and 160 mK curves at low bias- 
the voltage at which the zero bias dip starts developing - 
is approximately 0.6 fiY, much smaller than Ec/e, but in 
very good agreement with the effective Thouless energy 
E*c/e defined earlier. For nonlocal measurements on con- 
ventional NSN devices near zero bias, CAR is expected to 
give a negative contribution to the nonlocal resistance [2] , 
as two electrons with energies less than A, one from each 
normal metal, combine to form a Cooper pair. The de- 
crease in the nonlocal resistance in our samples suggests 
a similar process is happening in these samples, since the 
pair coupling in the proximity coupled normal metal is 
finite. Thus, two electrons, one from each normal lead, 
combine to form a correlated pair in the proximity cou- 
pled normal metal. Exactly how this process occurs is not 
clear, as it is not described by our current understanding 
of nonequilibrium transport in proximity-coupled normal 
metals. 

In summary, measurements on proximity-coupled nor- 
mal metals reveal a signature of long-range nonlocal 
quasiparticle correlations that may be related to the for- 
mation of pair correlations in the proximity-coupled nor- 
mal metal. Further study is required to elucidate the 
origin of these correlations. 
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PACS numbers: 



I. NONLOCAL MEASUREMENTS ON A NSN STRUCTURE AND CONTROL 
SAMPLES 



Figure Sla shows an image of another sample that was measured. The length L of the 
normal wire between the two superconductors is shorter than that of the sample discussed 
in the main text. This sample had additional leads to perform control measurements. Fig- 
ure Sib shows a nonlocal resistance measurement corresponding to earlier measurements on 
NSN structures.^ (the numbers refer to the contacts numbered in Fig. Sla, with the same 
convention to specify the four terminal resistance as in the main text). It can be seen that 
the shape of the curve is very similar to what was observed previously. In particular, as 
shown in the inset, there is no sharp dip in the nonlocal resistance near zero bias. To em- 
phasize the difference between this measurement in the NSN configuration and the nonlocal 
measurements on proximity-coupled normal metals, Fig. Sic shows the nonlocal resistance 
similar to Fig. 2 in the main text for the sample of Fig. Sla. A very pronounced dip is seen 
at zero bias. 

To demonstrate that the nonlocal resistance is observed only when there is a supercurrent 
between the two superconducting contacts, we show nonlocal resistance measurements on 
another sample that is nominally identical to the sample of Fig 1 of the main text, except 
that it does not have a second superconducting contact (Fig. S2a). It can be seen that the 
nonlocal resistance is zero, showing that the presence of two superconductors is required to 
see the effect. 

II. TEMPERATURE DEPENDENCE OF THE ZERO BIAS NONLOCAL RESIS- 
TANCE 

Figure S3 shows the temperature dependence of the three different zero bias nonlocal 
resistances of the sample of Fig. 1 in the main text. It should be noted that while this 
cooling temperature sweep was performed over 5 hours, there is still some lag between the 
temperature read by the mixing chamber thermometer and the actual temperature of the 
sample, in contrast to the data shown in Fig. 2a of the main paper. Nevertheless, these data 
show that the nonlocal resistance is essentially zero above 0.8 K and increases rapidly below 
this temperature. This increase is associated with the establishment of Josephson coupling 
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FIG. 1: (a) Scanning electron micrograph of another nonlocal sample, with terminals numbered 
as in the main text. The light areas are gold, the darker lines are Al. Size bar is 500 nm. (b) 
Nonlocal resistance measmement of the NSN configuration on the left, corresponding to previous 
measurements on NSN samples.^ The numbers refer to contacts used for current and voltage leads, 
as described in the main text. Inset: Expanded view of the zero bias region, showing that there 
is no dip in nonlocal resistance near zero bias, (c) Nonlocal resistance of the proximity-coupled 
normal metal for the sample of (a). Inset: Expanded view of the zero bias region, showing that 
there are pronounced dips in nonlocal resistances near zero bias. 



between the two superconductors, as discussed in the main text. The decrease in resistance 
at lower temperatures is consistent with the development of the dip in nonlocal resistance 
seen in the nonlocal dV/dl at low temperatures (Fig. 4 in the main text). 



III. NUMERICAL SIMULATIONS BASED ON THE QUASICLASSICAL THE- 
ORY OF SUPERCONDUCTIVITY 

The numerical simulations in this paper were based on the Usadel equations of the quasi- 
classical theory of superconductivity, using the Riccati parametrization. We used the public 
domain numerical solvers developed by Pauli Virtanen^, based on earlier work^"^. The 
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FIG. 2: (a) Scanning electron micrograph of control sample similar to that of Fig. 1 of the 
main text, but with no second superconductor. Size bar is 200 nm. (b) Nonlocal resistance 
measurement on the sample of a, showing that there is no nonlocal signal in the absence of the 
second superconductor. 
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FIG. 3: (a) Temperature dependence of the zero bias nonlocal resistance of the sample of Fig. 1 
of the main text. 

starting point for the simulations is the equation for the total current 



j(R,T) = eN^D / dE[M^^{d^hT) + QHl + MosI^r/i^)]. 



(1) 

Here Nq is the electronic density of states at the Fermi energy, D is the electronic diffusion 
coefficient, hx and hi are the tranverse and longitudinal quasiparticle distribution functions. 
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which in equihbrium at energy E in reservoirs at a potential V have the form 

(2) 



, 1 



, , E + eV\ , fE-eV 
tanh — — — ± 



2kBT J \ 2kBT 



the spectral supercurrent Q is given by 



g = 3? (2 V7 - 7 V7] ) (3) 
and the dimensionless diffusion coefficients M by 

M33=|iVp(|7r + l)(l7r + l) (4) 

and 

Mo3 = |iVp(|7|2-|7|2), (5) 

where = (1 + 77)"^. 7, 7 are the Riccati parametrization parameters that satisfy the 
coupled equations 

DV^-f-2N^\V-f\^ + 2iE^^0, (6) 



and 



L>V^7 - 2iV7| V7r + 2i£;7 (7) 



in the normal metal wires. The boundary conditions for the differential equations are that 

7 and 7 are zero at a normal reservoir, while on a superconducting reservoir 

A A* 
7« = - 7« = - (8) 

E + iy/\A\^ - {E + iS)^' E + iy/\A\^ - {E + i5)^' 

where A = |A|e"^ is the complex order parameter in the superconducting reservoir, being 
the phase of the superconductor. 

The first term on the RHS in Eqn (1) is the quasiparticlc current, the second term is the 
supercurrent, and the third term corresponds to conversion of quasiparticlc to supercurrent, 
and is typically negligible in a normal metal. At the nodes, the Riccati parameters are 
continuous, and their derivatives sum to zero. The spectral charge and energy currents also 
sum to zero along a node, and are conserved in a normal wire. 

To calculate the differential resistance as a function of the current, we use the sample 
geometry in Fig. 3a of the main text. A voltage V is applied to the normal contact Ni with 
the superconducting reservoirs 5*1 and 5*2 being at zero potential, and the resulting current 
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Iqp flowing from A^i is calculated under the boundary condition that no current flows into 
contacts N2, N3 and S2. In order to satisfy these boundary conditions, the voltages Vnii 
and Vni2 on the normal contacts N2 and A^3 and the phase difference (f) between the two 
superconducting contacts Si and S2 are adjusted in an iterative loop. As Iqp increases, a 
quasiparticle current Iqp2 flows into S2, which is counterbalanced by a supercurrent Is that 
flows from 5*2 to Si. As Is approaches the critical current, the simulations have greater 
difficulty converging, and we have shown in the figures in the main text only that range of 
Iqp over which the boundary conditions are satisfied. For the plot, the nonlocal differential 
resistances dVnii/dlqp and dVnn/dlqp are calculated numerically. Taking the length of the 
wire between the superconducting reservoirs to be L, the nonlocal resistances shown in the 
simulations of Fig. 3b of the main text correspond to the positions of the normal leads 
in the actual sample (Fig. la of the main text), with the length of each normal reservoir 
from the proximity coupled normal wire being 0.75 L. We have also performed simulations 
with different values of this length, and this value most closely resembles the shape of the 
experimental curve. The temperature of the simulations correspond to 20 mK. 

For the temperature dependence, a similar simulation was done, but only a small voltage 
±y was applied to A^i, and the resulting values used to calculate the differential resistances 
as a function of temperature. 
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